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Abstract:  In  eddy  current  non-destructive  testing  of  a  multi-layered  riveted  structure,  rotating  current 
excitation,  generated  by  orthogonal  coils,  is  advantageous  in  providing  sensitivity  to  defects  of  all 
orientations.  However,  when  used  with  linear  array  sensors,  the  exciting  magnetic  flux  density 
(Bx)  of  the  orthogonal  coils  is  not  uniform  over  the  sensor  region,  resulting  in  an  output  signal 
magnitude  that  depends  on  the  relative  location  of  the  defect  to  the  sensor  array.  In  this  paper, 
the  rotating  excitation  coil  is  optimized  to  achieve  a  uniform  Bx  field  in  the  sensor  array  area  and 
minimize  the  probe  size.  The  current  density  distribution  of  the  coil  is  optimized  using  the  polynomial 
approximation  method.  A  non-uniform  coil  design  is  derived  from  the  optimized  current  density 
distribution.  Simulation  results,  using  both  an  optimized  coil  and  a  conventional  coil,  are  generated 
using  the  finite  element  method  (FEM)  model.  The  signal  magnitude  for  an  optimized  coil  is  seen  to 
be  more  robust  with  respect  to  offset  of  defects  from  the  coil  center.  A  novel  multilayer  coil  structure, 
fabricated  on  a  multi-layer  printed  circuit  board,  is  used  to  build  the  optimized  coil.  A  prototype 
probe  with  the  optimized  coil  and  32  giant  magnetoresistive  (GMR)  sensors  is  built  and  tested  on 
a  two-layer  riveted  aluminum  sample.  Experimental  results  show  that  the  optimized  probe  has  better 
defect  detection  capability  compared  with  a  conventional  non-optimized  coil. 

Keywords:  eddy  current;  non-destructive  testing;  giant  magnetoresistive  sensor;  optimization 


1.  Introduction 

The  detection  of  deep,  embedded  cracks  in  multilayer  riveted  structures  [1-4]  is  a  major  challenge 
in  eddy-current  (EC)  non-destructive  testing  (NDT).  Ultra-low  frequency  excitation  along  with  giant 
magnetoresistive  (GMR)  sensor  [5]  to  measure  the  induced  magnetic  field  directly  (EC-GMR)  have 
been  used  to  increase  penetration  depth  as  well  as  guarantee  a  good  signal  to  noise  ratio  [6-8]. 

The  EC-GMR  technique,  using  a  single  linear  excitation  coil,  has  been  presented  to  detect  fatigue 
cracks  around  fasteners  in  multilayer  structures  [9-11],  where  the  induced  currents  are  primarily 
generated  along  one  direction.  A  linear  array  of  GMR  sensors  located  on  the  symmetry  plane  of  the 
coil  picks  up  the  normal  component  ( Bz )  of  a  magnetic  field  associated  with  eddy  currents.  This  design 
has  high  sensitivity  when  the  crack  is  perpendicular  to  the  current  direction.  However,  when  a  crack  is 
parallel  to  the  EC  flow,  the  small  disturbance  to  eddy  current  flow  reduces  the  detection  capability. 
In  order  to  address  this  problem,  an  excitation  method  that  generates  a  rotating  eddy  current  was 
presented  by  Yang  et  al.  [12].  The  rotating  current  excitation  is  produced  by  two  orthogonal  linear  coils 
oriented  in  x  and  y-direction,  as  shown  in  Figure  la.  Conventionally,  the  wires  of  the  coils  are  uniformly 
distributed,  which  is  referred  to  as  a  conventional  coil  in  the  following  text.  The  exciting  current  in  the 
x-direction  is  expressed  as  Ix  =  Iqco s  (cot) ,  while  the  current  in  y-direction  is  Iy  =  Iqc os  (cot  +  90°), 
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where  Jo  is  the  current  amplitude,  cu  angular  frequency.  The  induced  eddy  current  in  the  sample 
generated  by  the  orthogonal  coils  is  rotating  in  space  and  hence  sensitive  to  defects  in  any  orientation. 
A  probe  configuration  with  rotating  current  excitation  has  similar  sensitivity  to  cracks  emanating 
around  fastener  sites  in  all  orientations.  However,  with  two  orthogonal  coils,  there  is  only  one  point 
of  symmetry  at  the  center  of  the  coils.  The  normal  component  of  induced  flux  density  is  zero  at  the 
symmetry  point  in  the  absence  of  a  defect. 

When  using  a  probe  with  a  single  sensor  located  at  the  symmetry  point,  a  2D  raster  scan  is  needed 
that  makes  an  inspection  of  a  sample  surface  slow.  However,  if  an  array  of  sensors  is  used  to  increase 
inspection  speed,  the  sensor  array  will  be  located  on  the  symmetry  plane  of  one  coil  but  perpendicular 
to  the  symmetry  plane  of  the  other  coil.  It  leads  to  two  major  issues,  which  give  rise  to  some  challenges 
in  experimental  implementation,  namely,  (1)  a  strong  background  field  at  sensors  located  off  the 
center  point  and  (2)  non-uniform  induced  currents  in  the  test  sample  below  the  sensor  array  location. 
Even  theoretically,  these  issues  can  be  addressed  by  using  an  extremely  large  size  exciting  coil;  the  total 
size  of  the  probe  is  preferred  to  be  small  in  field  application.  Consequently,  it  is  necessary  to  make  the 
rotating  current  EC-GMR  probe  have  small  size,  as  well  as  robust  performance. 

The  differential  measurement  method  was  proposed  to  address  the  first  issue  of  eliminating  the 
strong  background  field  during  the  scan  in  [13].  In  this  paper,  the  excitation  coil  design  is  optimized  to 
generate  a  uniform  excitation  field  so  as  to  address  the  second  issue.  The  performance  of  the  optimized 
probe  design  is  analyzed  and  compared  with  a  conventional  probe. 


2.  Conventional  Uniform  Exciting  Coil 

As  shown  in  Figure  la,  in  the  rotating  current  configuration  using  orthogonal  coils,  the  linear 
sensor  array  is  located  on  the  symmetry  plane  of  the  x-direction  coil,  and  is  perpendicular  to  the 
symmetry  plane  of  the  y-direction  coil.  The  magnetic  flux  density  components,  Bx,  generated  by  the 
y-direction  coil,  is  not  uniform  over  the  sensor  region,  which  implies  that  the  eddy  current  generated 
by  the  excitation  coil  in  the  test  sample  is  not  uniform  over  the  sensor  region. 


GMR  array  y-direction  coil 


-►  x  x-direction  coil 


(a) 


x  (mm) 

(b) 


0  Unit  (mm) 

(c) 


Figure  1.  (a)  Top  view  of  orthogonal  coils  and  linear  array  sensors;  (b)  Simulation  result  of  response 
signals  due  to  defect  located  at  x  =  0  mm  and  x  =  30  mm;  (c)  Lateral  view  of  the  model. 


The  eddy  current  density  at  the  center  of  the  sensor  array  is  higher  in  magnitude  compared  to 
the  eddy  current  density  at  the  edge  of  the  sensor  array  at  the  same  depth.  This  implies  that  a  defect 
located  under  the  sensor  at  the  center  is  detected  with  a  stronger  response  compared  to  the  signal 
produced  at  the  sensor  when  the  defect  is  located  near  the  coil  edge  at  same  depth.  The  ID  scan  signals 
for  the  two  cases  are  shown  in  Figure  lb,  where  the  geometry  dimensions  of  the  model  are  as  shown 
in  Figure  lc.  When  the  probe  is  used  to  inspect  rows  of  rivets,  the  rivet  image  will  be  distorted  due  to 
this  inconsistency. 
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3.  Exciting  Coil  Optimization 

In  order  to  address  the  problems  described  above  and  improve  the  performance  of  the  probe, 
the  excitation  coil  design  is  optimized  to  make  the  excitation  field  uniform  across  the  sensor  array. 
The  lateral  view  of  the  y-direction  excitation  coil  and  the  sensor  array  is  shown  in  Figure  2.  The  coil 
consists  of  the  active  region  with  excitation  current  and  return  paths.  The  active  part  is  defined  by 
x  G  (— a ,  a)  while  the  return  paths  lie  in  the  interval  x  G  (b,  c )  U  (— c,  —b),  here  0  <  a  <  b  <  c. 


.  Sensor  array  A  z 

Air  gap  ...  Air  gap 


c  -b  -a 

:  x 

0 

a  b  c 

Return  path 

active  part 

Return  path 

Figure  2.  Lateral  view  of  the  y-direction  exciting  coil  and  sensor  array. 


Let  }(x)  (— c  <  x  <  c)  be  the  current  density  distribution  of  the  coil: 

J(x)  =  J(-x) 

[  J  (x)  dx  =  0 
Jo 


(1) 


(2) 


Here  Equation  (1)  is  due  to  the  symmetric  structure  of  the  coil.  The  current  density  distribution  is 
an  even  function  of  x  and  hence  it  is  necessary  to  optimize  only  for  (x  >  0).  Equation  (2)  indicates 
that  the  total  return  current  is  equal  to  the  total  excitation  current  for  each  half  of  the  coil  to  satisfy  the 
current  continuity  theorem. 

Since  the  return  part  is  relatively  far  away  from  the  sensor  array  and  the  width  of  the  return 
part  is  kept  as  small  as  possible  to  make  the  coil  size  small,  it  is  reasonable  to  assume  that  the  current 
distribution  in  the  return  path  is  constant  ( Jr ).  The  current  density  distribution  in  the  active  region  is 
a  function  of  x  and  needs  to  be  optimized.  Let  /(x)  be  written  as  nth  order  polynomial  in  the  region 
(x  G  [0 ,d\): 

YY  Pk  xk  0  <  x  <  a 

JW  =  <  o"°  a<x<b  <3> 

Jr  b  <  X  <  C 

Here  p (k  =  0, 1,2,  •  •  •  ,n)  is  the  factor  of  the  power  series.  Substituting  Equation  (3)  into 
Equation  (2): 


j  1  ^  * 

Jr**b=~c  ^Pkk- 

V  L  k=0  K 


k+1 


1 


(4) 


Since,  for  riveted  structure  inspection,  the  exciting  frequency  is  ultra-low,  e.g.,  100  Hz,  Bx  at  any 
observation  point  (x',  0,  z ')  in  air  is  written  as: 


Bx  =  j_G(x'  -x,z')]{x) 

where  G  (x,z)  is  Green's  function  expressed  as  [14,15]: 


dx 


(5) 


G  (x,z)  = 


yo  Z 


2 

I 


dy 


4nVx2  +  z2  J  x2  +  z2  +  y- 


(6) 


where  l  is  the  length  of  the  current  in  y-direction.  As  the  orthogonal  coil  is  square,  l  =  2c.  Invoking 
Equations  (l)-(5): 
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Bx  (x',z')  —  [G  (V  -  x,z')  +  G  (x'  +  x,z')]J  (x)  dx  (7) 

Substituting  from  Equations  (3)  and  (4)  in  Equation  (7)  and  interchanging  summation  and 
interpolation  we  have: 

Bx(xr, z'H  i 
k= 0 

We  define: 


vA  Cf  fG  fx'  -  x,z')  +  G  (x'  +  x,z')  lx*:rfx+  Cf 


r  /-?  /  / 


7k  (X z')  =  f* [G  (x'  -  x,  z')  +  G  (x'  +  x,  z')  Jx^dx  +  [G  (x'  -  x, z')  +  G  (x'  +  x, z')]  ^  ^  (9) 

such  that: 


Bx  ( x',z ')  =  £7^ 
k=0 

Suppose  there  are  m  observation  points  {x'i,  0,  z'/),  /  =  1,  2,  •••,  m.  Define  vector: 

Bx  =  [Bxi/  Bx2/  •  •  •  -^xm] 

Then  Equation  (10)  can  be  written  in  compact  form  as: 


(10) 


(11) 


Bx  =  YV 


(12) 


Here  7  is  a  (n  +  1)  x  m  matrix  that  is  defined  as: 

T  (i'j)  T i{x'j,zrj) 

P  is  the  unknown  power  series  factors  vector, 

P  =  [VO’Vl  ■■■Pn]T 

To  make  Bx  uniform  across  all  the  observation  points: 

7TP  =  V  =  [v,  v . . .  v\ T 


(13) 


(14) 


(15) 


where  v  is  a  non-zero  constant  value. 


3.1.  Size  of  Active  Region  "a" 

Equation  (15)  is  consistent  if,  and  only  if,  the  coefficient  matrix  7  and  the  augmented  matrix  have 
the  same  rank.  The  solvability  of  Equation  (15)  depends  on  a  number  of  design  parameters.  One  of 
the  most  important  parameters  is  the  size  of  the  active  region.  Consider  a  sensor  array  with  32  GMR 
sensors  located  at  xf  =  =b  (2/  —  1)  mm  (1  <  i  <  16).  Since  the  defects  are  typically  assumed  to  be  in 
the  second  layer  at  a  depth  of  6  mm  in  the  multiple-layer  structure;  it  is  required  to  have  a  uniform 
eddy  current  at  this  depth  to  avoid  the  distortion  due  to  non-uniform  excitation  fields.  The  observation 
points  where  the  excitation  field  is  measured  are  selected  as: 


i  =  1 

mm 

1  <  i  <  16 


(16) 


z'i  =  -6  1  <  i  <  17  mm  (17) 

As  there  are  17  observation  points  ( m  =  17),  n  is  set  to  be  16  to  make  7  a  square  matrix. 
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The  width  of  the  return  path  should  be  one  quarter  of  the  active  region  to  guarantee  that  there  is 
enough  space  for  return  wires  when  the  coil  is  fabricated  on  a  printed  circuit  board  (PCB). 

c~b  =  \  (18) 

Assume  the  air  gap  width  between  the  active  part  and  the  return  path  (wajr)  is  5  mm.  Then 

b  —  a  =  5  mm  (19) 


The  effect  of  "a"  is  studied  numerically  by  sweeping  parameter  " a "  from  25  mm  to  50  mm  at 
increments  of  1  mm.  For  each  value  of  “a” ,  “b”  and  “c”  are  calculated  according  to  Equations  (18) 
and  (19),  7  is  calculated  from  Equations  (9)  and  (15)  is  analyzed.  It  is  found  that  Equation  (15)  is  not 
consistent  when  a  <26  mm.  When  a  >27  mm,  a  non-zero  solution  of  P  is  derived  from  Equation  (15) 
from  which  /  (x)  is  calculated  using  Equation  (3).  Some  of  these  plots  of  /  (x)  are  presented  in  Figure  3a. 
It  is  seen  that  the  plot  of  /  (x)  oscillates  between  positive  and  negative  values  when  a  <  32  mm.  So  "a" 
should  be  no  less  than  33  mm,  which  is  reasonable  considering  that  the  sensor  array  is  located  at 
x  G  (—31,  31).  Figure  3b  presents  the  optimized  current  density  distribution  when  33  <  a  <  39, 
in  which  case  /  (x)  is  always  positive. 


Figure  3.  Optimization  result:  current  density  distribution  vs.  x  for  different  values  of  parameter  "a". 
(a)  27<  a  <  45;  (b)  33  <  a  <  39,  which  is  the  possible  range  of  " a ”  for  optimized  probe. 


Smaller  values  of  “a”  are  preferred,  since  the  overall  probe  size  is  preferred  to  be  small.  However, 
it  should  be  noticed  that  when  “a”  is  small,  the  variation  of  /  (x)  is  large,  e.g.,  when  a  -  33  mm, 
the  maximum  current  density  is  3248  A/m,  while  the  minimum  current  density  is  3.78  A/ m.  Such  a 
distribution  of  /  (x)  cannot  be  realized  on  a  PCB.  So  a  constraint  is  introduced  as  follows: 


max  |/  (x)  | 
min  |7  (x)  | 


(20) 


Keeping  fixed  wajr/  the  variation  of  rj  vs.  " a "  is  plotted  in  Figure  4. 
we  need  a  >  36  mm. 


For  rj  to  be  less  than  10, 


Figure  4.  Calculation  results:  tj  vs.  parameter  "a"  for  fixed  air  gap. 
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3.2.  Effect  of  Air  Gap  Width  "wair" 

The  effect  of  the  width  of  the  air  gap  between  the  active  region  and  the  return  path  (wajr)  is  studied 
in  this  section.  zvajr  is  varied  from  1  mm  to  9  mm  at  increments  of  1  mm.  The  value  of  p  defined  in 
Equation  (20)  is  calculated  for  different  values  of  wajr.  The  plots  of  rj  vs.  “a”  for  different  values  of  wajr 
are  presented  in  Figure  5. 

For  each  value  of  wa\r>  there  is  a  minimum  possible  “a”  (amin)  to  satisfy  p  <  10,  as  listed  in  Table  1. 
The  total  coil  size  (w)  can  be  calculated  according  Equation  (21). 

w  =  2.5a  +  2wajr  (21) 

As  the  air  gap  increases,  the  value  of  "amm"  decreases,  hence  a  good  tradeoff  is  provided  by 
allowing  air  gap  =  5  mm,  and  a  =  36  mm. 


Figure  5.  Calculation  results:  p  vs.  parameter  “a”  for  different  air  gap  width  wair 


Table  1.  Minimum  coil  size  satisfying  p  <  10  for  different  air  gap  width. 


Air  Gap  (mm) 

1 

2 

3 

4 

5 

6 

7 

8 

9 

"amin"  (mm) 

40 

39 

38 

37 

36 

36 

35 

35 

35 

w  (mm) 

102 

101.5 

101 

100.5 

100 

102 

101.5 

103.5 

105.5 

3.3.  Coil  Design 

Using  wajr  =  5  mm,  a  =  36  mm,  the  width  of  each  line  of  the  coil  dw  is  derived. 


dw  (x) 


1  x  max  !/(*)! 
/  0)  dm 


(22) 


where  dm  is  the  minimum  wire  width,  that  is  decided  by  the  minimal  wire  clearance  distance  of  PCB 
fabrication  processing  and  the  required  copper  wire  width  of  carrying  the  exciting  current.  Here  we  use 
dm  =  0.25  mm.  The  positions  of  each  line/wire  xWi  are  derived  according  to  the  following  equations. 


%ZVi  — 


dw  (0) 

2 


(23) 


xm+i  xwi 

The  final  design  of  optimized 
presented  in  Figure  6. 


dw  (fw,+1f  —  ^  ,  i>  1  and  xWj+1  <  a  (24) 

exciting  coil  with  50  turns,  non-uniform  distributed,  wires  is 
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Figure  6.  Optimized  coil  design  50  turns  non-uniform  distributed  wires:  (a)  wire  separation  (dw  )  vs.  x 
and  (b)  top  view  of  the  optimized  coil. 


4.  Model  Based  Performance  Analysis 

A  3D  finite  element  model  (FEM)  is  used  to  analyze  the  performance  of  the  optimized  coil  and 
compare  the  results  with  those  obtained  using  conventional  coil.  The  formulation  is  based  on  reduced 
magnetic  vector  potential  formulation  [16-18],  which  avoids  re-meshing  of  coil  for  each  position  during 
scanning.  Commercial  finite  elements  simulation  software  COMSOL  Multiphysics  (COMSOL,  Inc., 
Palo  Alto,  CA,  USA)  with  Matlab  LiveLink  is  used  to  build  the  model  and  do  numerical  calculation. 

The  magnetic  fields  in  air  generated  by  the  optimized  coil  and  conventional  coil  are  calculated. 
The  two  coils  have  identical  size  and  number  of  turns.  A  current  of  frequency  100  Hz  is  used  as 
excitation  input.  The  frequency  of  100  Hz  is  chosen  to  have  adequate  skin  depth  to  penetrate  the 
sample  thickness.  It  is  seen  from  the  results  presented  in  Figure  7 a  that  the  variation  of  Bx  in  the 
sensor  region  is  greatly  reduced  in  the  case  of  the  optimized  coil  in  comparison  to  that  obtained  using 
a  conventional  coil.  Although  the  normal  component  magnetic  field  ( Bz )  is  not  an  objective  function 
during  this  optimization,  as  seen  in  Figure  7b,  the  distribution  of  the  Bz  also  becomes  relatively  flatter 
with  the  optimized  coil,  in  comparison  to  that  obtained  with  a  uniform  coil. 


x  10'4 


(a) 


(b) 


Figure  7.  Numerical  simulation  results  of  (a)  Bx  at  the  observation  points  and  (b)  Bz  at  the  sensor  array 
for  the  cases  of  optimized  coil  and  uniform  coil. 


4.1.  Detection  of  Subsurface  Defects  at  Different  Locations 

Subsurface  defects  located  at  the  center  of  the  coil  (x  =  0  mm)  and  located  off  center  (x  =  30  mm) 
are  simulated.  The  sample  and  defect  geometry  dimensions  are  specified  in  Figure  8a.  Both  optimized 
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coil  and  conventional  coil  are  simulated.  A  normal  component  magnetic  field  ( Bz )  1  mm  above  the 
coil  is  measured.  To  highlight  the  defect  signal,  a  background  signal  is  obtained  by  placing  the  coil 
on  top  of  a  defect  free  aluminum  sample  and  subtracted  from  the  measured  signals.  The  simulation 
results  are  presented  in  Figure  8b.  For  conventional  uniform  distribution  coil,  the  peak  signal  is 
reduced  to  39.23%  when  the  defect  is  moved  from  the  center  of  the  coil  to  an  off-center  location  at 
x  =  30  mm  demonstrating  that  a  defect  located  near  the  coil  edge  is  detected  with  a  much  weaker 
response  compared  to  the  signal  produced  by  the  same  defect  at  the  coil  center.  However,  with 
the  optimized  coil,  the  peak  signal  only  reduced  to  81.25%.  In  the  model  based  study,  the  exciting 
current  was  kept  constant,  in  which  case  the  peak  signal  of  the  defect  located  under  the  center  of  the 
optimized  coil  is  smaller  than  that  under  the  center  of  the  conventional  coil;  but  it  doesn't  impact 
defect  detectability  of  the  sensor  since  the  signal  amplitude  can  be  easily  increased  by  using  a  larger 
exciting  current.  In  practice,  exciting  current  amplitude  is  tuned  to  make  the  strongest  signal  almost 
equal  to  the  measurement  range  to  obtain  best  sensitivity.  In  summary,  it  is  seen  from  the  simulation 
result  that  the  optimized  coil  has  better  performance  for  detecting  a  defect  located  under  the  edge  of 
the  sensor  array. 


Return  part  Exciting  part  _  Return  part 


2  Layers  aluminum  sample 

Defect  2 x2>^  5  30  3641  50  (mm) 


(a) 

x  IQ'7 


(b) 


Figure  8.  (a)  Lateral  view  of  tested  sample  and  excitation  coil;  (b)  Simulation  results  of  the  optimized  coil 
and  uniform  distribution  coil  for  a  defect  located  at  x  =  0  mm  and  x  =  30  mm  after  background  subtraction. 

4.2.  Detection  of  Defect  under  Fastener  Head 

Next,  defect  and  fastener  geometry  is  modeled  using  the  finite  element  model.  The  test  specimen 
is  a  two-layered  aluminum  structure  as  shown  in  Figure  9.  A  radial  notch  spans  the  thickness  of 
the  second  layer.  Rotating  current  excitation  with  two  identical  orthogonal  coils  is  used  in  these 
simulations.  The  excitation  frequency  is  100  Hz  and  current  amplitude  is  1  A.  The  orthogonal  coils 
are  located  at  a  1  mm  liftoff  from  the  top  surface  of  the  test  sample.  A  C-scan  image  is  produced  by 
calculating  the  normal  component  of  the  magnetic  field  ( Bz )  values  at  a  1  mm  liftoff  distance  from  the 
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coils  at  each  scan  point.  To  highlight  the  rivet  and  defect  signal,  the  background  signal  obtained  by 
placing  the  coils  on  top  of  defect  free  aluminum  sample  is  subtracted  from  the  measured  signal. 


z 


Figure  9.  Geometry  of  the  test  sample  used  in  the  finite  element  model  (FEM)  model. 


Two  different  cases  are  simulated:  (1)  the  rivet  center  is  aligned  with  the  orthogonal  coils' 
center  (x  =  0);  (2)  the  rivet  center  is  shifted  by  20  mm  along  x  axis  (x  =  20).  The  simulation  results 
showing  magnitude  of  the  normal  component  magnetic  field  ( Bz )  is  calculated  and  plotted  in  Figure  10. 
The  improved  performance  of  the  optimized  coil  is  seen  visually  in  the  C-scan  image.  The  improvement 
in  performance  can  be  quantitatively  estimated  as  follows. 


Figure  10.  Simulation  results  for  cases:  (a)  conventional  coil,  rivet  center  located  at  x  =  0  mm; 
(b)  conventional  coil,  rivet  center  located  at  x  =  20  mm;  (c)  optimized  coil,  rivet  center  located  at 
x  =  0  mm;  (d)  optimized  coil,  rivet  center  located  at  x  =  20  mm. 


To  quantify  the  distortion  introduced  by  the  presence  of  a  defect,  all  contiguous  pixels  above 
a  preselected  threshold  value  v  are  identified  and  a  contour  plot  enclosing  the  region  containing  pixels 
is  drawn.  The  contour  plot  obtained  by  setting  v  equal  to  30%  of  the  peak  value  is  plotted  in  Figure  10. 
Defining  the  deformation  coefficient  e  as  [13]: 


rmax  r0 

ro 


x  100% 


(25) 


where  rg  is  the  radius  of  the  contour  plot  of  the  image  obtained  with  a  defect  free  rivet.  rmax  is  the 
maximum  distance  from  the  rivet  center  to  the  contour  plot.  The  deformation  coefficient  is  a  simple 
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straightforward  parameter  for  evaluation  of  the  distortion  of  the  signal  in  a  rivet  area.  Some  other  more 
complex  algorithms  may  be  applied.  In  fact,  our  lab  colleagues  tried  different  methods.  For  example, 
J.  Kim  et  al.  used  principal  components  analysis  to  classify  the  data  [19],  a  multi-frequency  mixing 
algorithm  was  utilized  for  the  suppression  of  the  rivet  signal  [20].  As  this  paper  is  mainly  about  sensor 
optimization,  the  deformation  coefficient  is  used  to  quantize  the  signal.  The  deformation  coefficients 
of  the  simulation  results  are  compared  in  Figure  11. 

It  is  found  that  when  the  rivet  center  is  shifted  by  20  mm  from  the  coil  center  along  x  axis,  the 
distortion  coefficient  for  uniform  distributed  coil  drops  from  32.42%  to  12.41%.  However,  e  for  the 
optimized  coil  only  reduces  from  34.8%  to  20.8%.  The  optimized  coil  has  better  defect  indication  when 
the  rivet  and  defect  are  misaligned  with  the  coil  center. 


40 


30 

g 

|  20 

C/3 
Q_ 

LU 

10 

0 

Figure  11.  Simulation  Results — Comparison  of  reduction  in  e  of  the  conventional  coil  and  optimized 
coil  for  centered  and  off-center  rivets. 

5.  Experimental  Validation 

5.1.  Multi-Layer  Coil  Design 

A  prototype  probe  with  the  optimized  coil  is  fabricated  and  instrumented  to  test  the  performance 
of  the  proposed  design.  The  multi-layer  coil  is  fabricated  using  a  multi-layers  print  circuit  board 
(PCB).  A  novel  coil  design,  named  multi-layer  coils  (MLC),  is  proposed  to  fabricate  the  optimized  coil. 
Each  coil  consists  of  at  least  two  layers.  Coppers  from  different  layers  stagger  to  cover  the  air  gaps, 
as  shown  in  Figure  12,  so  that  the  current  distribution  is  more  continuous  in  comparison  with  that  of 
a  single  layer  coil. 


x=0mm  x=20mm 


0.3mm 


Figure  12.  Schematic  of  the  lateral  view  of  the  multi-layer  coil  design. 


5.2.  Experiment  Setup 

The  experimental  system  is  shown  schematically  in  Figure  13a.  GMR  Sensor  Model  GF708  from 
Sensitec  (Lahnau-Waldgirmes,  Germany)  is  used  to  build  the  prototype.  A  linear  sensor  array  with 
32  sensors  is  used.  The  distance  between  two  sensors  is  1.6  mm.  The  normal  component  of  the 
magnetic  flux  density  Bz  is  measured  by  the  array  sensors.  The  sensor  array  is  located  on  the  symmetry 
plane  of  the  x-direction  coil,  1  mm  above  the  coil. 

A  sinusoidal  current  of  amplitude  150  mA  and  frequency  100  Hz  is  applied  to  the  two  orthogonal 
coils  with  a  90°  phase  difference.  A  lock-in  amplifier  model  RF840  from  Stanford  Research  System  is 
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used  to  obtain  the  baseband  signal  with  the  best  possible  signal  to  noise  ratio  (SNR).  The  output  of  the 
lock-in  amplifier  is  sampled  and  stored.  The  experimental  measurements  are  obtained  using  a  high 
resolution  scanning  system. 


Orthogonal  coils  and 
double  sensor  arrays 


1 

Multiplexer  and 

1 

preamplifier 

—  Filter  -  Amplifier  -  Lock-in 


X  X  Power  amplifier 

I — - \9(f\ - 


DDS 


ADC 


Scanner 


Figure  13.  Block  diagram  of  experiment  system. 


5.3.  Experimental  Results  and  Discussion 

The  performance  of  the  optimized  probe  is  compared  with  that  of  a  conventional  coil  for  the  cases 
of  the  rivet  located  at  x  =  20  mm  off  the  coils'  center.  A  defect  free  aluminum  fastener  and  a  fastener 
with  a  6  mm  length  defect  in  the  second  layer  are  scanned.  The  experimental  results  are  presented 
in  Figure  14,  in  which  the  contour  selected  at  30%  of  the  maximum  value  in  the  image  is  plotted. 
Theoretically,  the  image  of  the  defect  free  rivet  should  be  circular  [12],  as  seen  in  Figure  14a.  However, 
it  is  seen  from  Figure  14b  that  the  image  of  the  defect  free  rivet  generated  using  the  conventional  coil  is 
distorted  due  to  the  non-uniform  distribution  of  the  excitation  field  and  eddy  currents. 


Figure  14.  Experiment  results  of  (a)  defect  free  rivet  generated  using  optimized  coil;  (b)  defect  free 
rivet  generated  using  conventional  coil;  (c)  rivet  with  a  6  mm  length  2nd  layer  defect  generated  using 
optimized  coil;  (d)  rivet  with  a  6  mm  length  defect  generated  using  conventional  coil. 
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The  deformation  coefficient  e  defined  in  Equation  (25)  is  calculated  to  quantify  the  performance, 
as  shown  in  Figure  15.  For  optimized  coil,  e  is  6.75%  for  the  defect  free  rivet  and  increases  to  18.8% 
for  the  rivet  with  a  6  mm  length  defect.  However,  for  the  conventional  coil,  e  changes  from  21.8%  to 
23.6%  for  the  defect  free  and  defective  rivet.  The  optimized  coil  has  better  defect  detection  capability 
compared  with  the  conventional  coil. 


Non-defect  6mm  defect 


Figure  15.  Experimental  Results — Comparison  of  reduction  in  e  of  the  conventional  coil  and  optimized 
coil  for  centered  and  off-center  rivets. 

6.  Conclusions 

In  this  paper,  the  current  distribution  of  rotating  current  excitation  is  optimized  with  respect 
to  a  linear  sensor  array.  A  non-uniform  coil  design  was  derived  to  generate  the  optimized  current 
distribution.  The  main  desirable  feature  of  the  optimized  coil  design  is  that  the  excitation  field  is  more 
uniform  in  the  sensor  array  region  than  that  of  the  conventional  coil. 

The  performance  of  the  optimized  design  was  studied  numerically  based  on  FEM  models  and 
validated  experimentally  using  a  prototype  probe  to  test  a  riveted  multi-layer  aluminum  sample. 
Simulation  results  show  that  for  the  conventional  coil,  a  subsurface  defect  located  near  the  coil  edge  is 
detected  with  a  much  weaker  response  compared  to  the  signal  produced  when  the  defect  is  located  at 
the  coil  center.  However,  the  reduction  in  signal  magnitude  is  much  less  in  the  case  of  an  optimized  coil. 
For  the  inspection  of  a  defect  located  under  the  rivet  head,  the  distortion  coefficient  of  the  optimized 
coil  is  greater  than  that  of  the  uniform  coil  when  the  fastener  and  defect  are  not  aligned  with  the 
coils'  center. 

A  new  multi-layer  coil  structure  was  proposed  to  fabricate  the  optimized  coil  of  the  prototype 
probe.  Experiment  results  demonstrate  that  the  optimized  probe  has  a  higher  detection  capability  of 
a  subsurface  defect  in  a  multi-layer  structure  compared  to  a  conventional  coil. 
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